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Abstract 

We calculate the probability distribution of a single component of peculiar velocities 
due to cosmic strings, smoothed over regions with a radius of several Mpc. The 
probability distribution is shown to be Gaussian to good accuracy, in agreement with 
the distribution of peculiar velocities deduced from the 1.9 Jy IRAS redshift survey. 
Using the normalization of parameters of the cosmic string model from CMB mea- 
surements, we show that the rms values for peculiar velocities inferred from IRAS 
are consistent with the cosmic string model provided that long strings have some 
small-scale structure. 
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1 Introduction 



Cosmic strings might be responsible for the formation of large-scale structure which 
is observed in the Universe today. They are topological defects which form in a 
phase transition in the early universe. Topological defect models of structure forma- 
tion are an alternative to inflationary models. In inflationary models one assumes 
the existence of a scalar field which drives a phase of very rapid expansion in the 
early universe, and quantum fiuctuations produced during this phase turn into clas- 
sical density perturbations which grow by gravitational instability into present-day 
structures. In topological defect models the defects, which are present at all times 
including today (and could possibly be seen directly some day), act as seeds around 
which structures form. It is important to study theories from both classes of models 
in order to see which predictions can distinguish between different models, and which 
predictions are too generic, arising in models based on widely different assumptions. 

Among the topological defect models, cosmic strings were the first to be investi- 
gated in more detail (Brandenberger 1991), recently textures have also been subjected 
to closer scrutiny (Turok 1991). The cosmic string model looks promising so far, since 
normalizations of the parameters occuring in the model obtained from observation 
agree with each other (Bennet, Bouchct & Stcbbins 1992, Perivolaropoulos 1993a). 
These obervations include galaxy redshift surveys and measurements of the temper- 
ature anisotropics in the cosmic microwave background radiation (CMB). Moreover, 
cosmic strings naturally produce filaments and planar structures in the matter distri- 
bution (Silk & Vilenkin 1984, Vachaspati 1986, Stebbins et al. 1987), in encouraging 
agreement with recent galaxy redshift surveys. Another nice feature of the cosmic 
string model is that it works well in the context of hot dark matter (Brandenberger 
1987), for which there are candidates known to exist (the /i and r neutrinos), al- 
though they may be massless. In contrast, none of the candidates for cold dark 
matter required in infiationary models has been detected. 

Topological defect models have not been as closely investigated as inflationary 
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models, and more work is necessary to quantify the predictions of the cosmic string 
model. It is especially important to find features accessible to observations in which 
cosmic strings differ from the class of inflationary models, and statistics sensitive to 
these features. Recently, the probability distribution of peculiar velocities has been 
put forward as a statistic which might be capable of discriminating between cosmic 
strings and inflation (Kofman et al. 1994). Here we show that according to the 
calculations in a simple string toy model, this statistic cannot be used to differentiate 
between cosmic strings and inflation. 

The IRAS 1.9 Jy redshift survey was analyzed (Strauss et al. 1990 & 1992, 
Yahil 1991) to obtain a uniform galaxy-density map. The peculiar velocities were 
reconstructed from it via a self-consistent iterative scheme assuming linear biasing 
between the density fluctuations of galaxies and mass (Nusser et al. 1991). Using 
these results, Kofman et al. (1994) determined the probability distribution of a 
single component of peculiar velocities of regions smoothed over several Mpc, 
with the result that it is consistent with the underlying distribution p{vx) being a 
normal distribution. Their result for the Gaussianity of the velocity distribution is 
still tentative because of the hmited volume sampled, and because velocities were not 
measured directly but deduced from a redshift survey. 

Inflationary models predict a normal distribution for p{vx) , whereas some devia- 
tion from gaussianity is expected in the cosmic string model since individual strings 
impart coherent velocity perturbations over extended regions, as we will see in the 
following section. The question is, however, just how large this deviation from gaus- 
sianity is and whether it is big enough to be detected by current observations. For 
a certain region receives velocity perturbations from many strings, and by the cen- 
tral limit theorem many nongaussian perturbations can add up to a gaussian sig- 
nal. Scherrer (1992) has shown that for seed models the velocity field can be very 
nearly Gaussian even if the density field is nongaussian. In a model where randomly 
distributed point masses, all of the same mass, accrete matter gravitationally in a 
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universe dominated by hot dark matter, he found that a very low seed density of less 
than 10~^Mpc~^ is required in order for p{vx) to be nongaussian. 

The main aim of this paper is to calculate p{vx) within the cosmic string model in 
order to quantify its departure from gaussianity, and to find out if the IRAS results 
are also in agreement with this stringy probability distribution p{vx)- Previously, 3 
dimensional rms velocities have been obtained analytically within the cosmic string 
model (Vachaspati 1992, Perivolaropoulos & Vachaspati 1994), and numerical sim- 
ulations have been performed (Vachaspati 1992, Hara, Mahonen & Miyoshi 1993) 
to obtain the probability distribution of peculiar velocities, and it was found that 
the distribution is quite gaussian. Here we present the first analytical calculation of 
p{vx) , within a model for the string network which has been previously employed 
to make predictions about temperature anisotropies in the CMB (Perivolaropoulos 
1993a & 1993b, Moessncr, Perivolaropoulos & Brandenberger 1994) and the magni- 
tude of peculiar velocities (Vachaspati 1992, Perivolaropoulos & Vachaspati 1994). 
The dependence of the deviation from gaussianity on one of the parameters of the 
model, namely the number i/ of strings per Hubble volume in the strings' scaling 
solution, is shown explicitly. 

Our conclusion is that on scales of several Mpc, p{vx) from strings deviates 
only slightly from a normal distribution. On the smallest scales and for the smallest 
number of strings per Hubble volume, the largest deviation from gaussianity is ex- 
pected. Performing a -test of the stringy p{vx) for the smallest scale of 6/i~^ Mpc 
given in Kofman et al. (1994), with imaginary data binned in the same way as in that 
paper, but drawn from an underlying gaussian distribution (worst case), yields the 
result that the stringy probability distribution is in agreement with these data. This 
agreement holds for all values of v , the number of long strings per Hubble volume 
in the strings' scaling solution, including u — 1. This shows that the hope expressed 
in Kofman et al. (1994) that the scenario for the formation of large scale structure, 
where widely separated strings accrete matter in wakes behind them, can be ruled 
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out using the statistic p{vx) is not realized. 

In the next section we will describe the cosmic string model and the mechanism 
for the production of velocity perturbations. In the third section we will describe 
the analytic model for the string network and calculate the moment generating func- 
tion for a single component of the peculiar velocities averaged over several Mpc. 
Then we will obtain the probability distribution p{vx) from this moment generating 
function, and compare it with observations. 

2 Cosmic Strings and Structure Formation 

Cosmic strings are linear topological defects formed at a phase transition in the very 
early universe (Vilenkin 1985). Those originating in the symmetry braking of a Grand 
Unified Theory possess an enormous mass per unit length fi {Gfi ~ 10^^ , where 
G is Newton's constant), and they can be responsible for the structures observed 
today. Strings can have no ends, so they are formed as either infinitely long strings 
or closed loops. After formation the network of strings quickly evolves towards a 
scaling solution where the energy density in long strings remains a constant fraction 
of the total background energy density. This is achieved by intercommutations and 
self-intersections of the strings leading to the production of small loops, which then 
decay by emitting gravitational radiation. In this way, some of the energy input into 
the string network coming from the stretching of the strings due to the expansion of 
the universe is transferred to the background. Long strings are straight over distances 
of the order of the horizon, so that the scaling solution can be pictured as having a 
fixed number v of long strings per Hubble volume at any given time. 

Initially loops were thought to make the dominant contribution to structure for- 
mation. At distances much larger than their radius they act as point masses and ac- 
crete surrounding matter (Turok & Brandenberger 1986, Stebbins 1986, Sato 1986). 
Improved cosmic string evolution simulations (Bennet & Bouchet 1988, Albrecht & 
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Turok 1989, Allen & Shellard 1990) showed that more of the energy density is in 
long strings, and they are therefore more important, accreting matter in the form of 
wakes behind them as they move through space (Brandcnbcrger, Perivolaropoulos & 
Stebbins 1990, Perivolaropoulos, Brandenberger & Stebbins 1990). Spacetime around 
a long straight cosmic string can be pictured as locally flat, but with a deficit angle 
of SnGfj, (Vilenkin 1985). Therefore a string moving relativistically with velocity Vg 
imparts velocity perturbations to surrounding matter towards the plane swept out 
by the string. If small-scale structure is present on the string, there is in addition a 
Newtonian force towards the string. The magnitude of this velocity perturbation is 
given by (Vachaspati & Vilenkin 1991, VoUick 1992) 

u = 4nGi,^sVsf, / = 1 + (1) 

In the absence of small-scale structure on the string,the tension T of the string is 
equal to its mass per unit length /i, and f — 1- If small-scale structure is present, 
II denotes the mass per unit length obtained after averaging over the small scale 
structure, T ^ fx, and / > 1. We consider the perturbations caused by strings after 
teq, the time of equal matter and radiation, in a universe filled with hot or cold dark 
matter. By the present time, the initial velocity perturbation imparted to the dark 
matter at time ti has grown to (Brandenberger 1987, Stebbins 1987, Kara & Miyoshi 
1990) 

Ui ^ OAu^/^ (2) 

Due to compensation (Traschen, Turok & Brandenberger 1986, Veeraraghavan & 
Stebbins 1992), the deficit angle of strings which are straight over a horizon distance, 
extends out only to a distance of one Hubble radius from the string, so that 
matter which is farther away does not receive any velocity perturbations. The velocity 
perturbation given in eq.(l) is independent of distance from the string (up to the 
Hubble radius), so that cosmic strings impart coherent perturbations over regions of 
the size of half a Hubble volume. 
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3 Moment Generating Function 

The moment generating function (mgf) Mx{t) of a random variable X is defined by 

Mx{t) = (exp*^) (3) 

where the brackets denote the ensemble average, and it contains complete information 
about X. In the following, X — denotes the random variable for the component 
of the peculiar velocities smoothed over a region V of comoving radius in a fixed 
direction Cx ■ We will calculate the mgf of X in order to obtain the moments (and 
cumulants) and probability distribution p{vx) of X from it. The moments of X, 
II j —< X^ > are given by 

and the cumulants Cj are defined by 

The probability distribution p{vx) can be expanded in an asymptotic series called 
Edgeworth series (Scherrer & Bertschinger 1991, Stuart & Ord 1987) in terms of the 
quantities 

A, = cj/d^' (6) 
and Hermite polynomials Hn{x) defined by 

H4x) = (-l)«e^^/^£^e-V2 (7) 

For distributions with vanishing odd moments as in our case, the expansion is 
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p(5) = ^-^[1 + —HJ5) + —He(5) + '-^^^^^Hs(5) + • • •] (8) 
V27r ^ 24 ^ ^ 720 ^ ^ 40320 ^ ^ ^ ^ ^ 

where 5 — v^jo , and a is the standard deviation of X. Since it is an asymptotic 
expansion, the remainder is of the order of the last term included (Erdelyi 1956). 
The mgf has an important property which makes it useful for calculations. For 
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independent random variables X and Y, the mgf of the sum is the product of the 
individual moment generating functions 

MxMt) = Mx{t)MY(t) (9) 

The following calculations are carried out within an analytical model for the string 
network which has previously been used to obtain the temperature anisotropics in the 
CMB and the magnitude of peculiar velocities from strings (see references given in 
the introduction). According to the scaling solution for cosmic strings, there is a 
fixed number of long strings present per Hubble volume at any given time. After 
about one expansion time of the universe (Hubble time), t 2t, these strings will 
typically have self-intersected or intercommuted, so that the resulting strings are 
uncorrelated with the ones at the previous Hubble time. We will assume that during 
one expansion time u long strings move across the Hubble volume. Each string is 
assumed to be straight over one horizon distance, and the effect of all strings is 
taken to be the superposition of the effects of the individual strings. The fact that 
small-scale structure varies along strings is neglected, so that we might somewhat 
underestimate the degree of non-Gaussianness of the velocity distribution. We also 
assume that the strings' positions, velocities and orientations at each Hubble time are 
random and uncorrelated, although this is not strictly true, since - to mention one 
reason - the string network has the form of a self-avoiding random walk. According 
to the two previous assumptions, the random variable for the total peculiar velocity, 
X, is the sum of independent random variables for the velocity perturbations from 
the individual strings. By eq.(9) we can therefore reduce the calculation of the the 
mgf for X to that of the mgf for the contribution of only one string, and take the 
products afterwards. In fact we know that there arc u strings per Hubble volume 
on average. The products in eq.(9) simplify if we take a Poisson distribution for the 
number of strings per Hubble volume instead of assuming the presence of exactly u of 
them (Scherrer & Bertschinger 1993). We can picture this as having a reservoir of n 
strings, each with a probability u/n oi being present in a particular Hubble volume, 
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in the limit that n — > oo. So if My.(i) denotes the mgf for a single component of 
peculiar velocities of the region V due to one string present at time U in the region's 
Hubble volume, then 

MxAt) = lim[^(exp*^0 + (l--)r 

= exph(My,(i)-l)] (10) 

where Xi is the random variable for the contribution of all strings at Hubble time 
ti to the velocity perturbation of V, and Vi denotes the average number of strings 
having an effect on V at time ij, i.e. those strings which are within a distance of one 
Hubble radius of V at time U (see eq.(18)). 

Since X = J2iLi , where is the number of expansion times since t^q 

TV^log,^ (11) 

''eg 

and the X^ are assumed to be independent, 

Mxit) = flMxXt) 

i=l 

N 

= expEi.,(My,(t)-l)] (12) 
1=1 

We will now calculate the mgf for Yi, the random variable for the component in 
the fixed direction of the peculiar velocities of a region V of comoving radius R 
(the smoothing radius) due to one string affecting the region at time t,. We have 
to perform the ensemble average over positions , orientations and directions of 
velocity Vg of the string. 

For a long straight string only transverse velocities are observable, and we assumed 
is and Vg to be random unit vectors. Therefore the unit normal e — ig x Vs oi the 
plane swept out by the string, i.e. the direction in which matter receives velocity 
perturbations, is itself a random unit vector. Consequently, the projection s = e ■ Cx 
is uniformly distributed over the interval [—1, 1], and the magnitude of the velocity 



8 



perturbation from one string has to be multiplied by s to get the component in 
direction e-x- 

During one expansion time a string sweeps out a plane towards which matter 
within a distance of one Hubble radius receives velocity perturbations, which have 
grown to Ui (see eq.(2)) by today. The possible values for the random variable Fj, 
the projection of the peculiar velocity of V in direction Cx due to one string at time 
ti, are a function of the perpendicular distance r of the centre of V to this plane: 

Uiif) — sUir/Ri for < r < Ri 

yi{r) = sui for Ri < r < Hf^ - Ri 

y^{r) = su, "''~^^~^'^ for Hr^-R,<r< Hr^ + R, (13) 

Ri is the physical size of the comoving radius R at time ti. Strings within a distance 
of H~^{ti) = of the region V can affect it. We distribute the centres C of these 
planes randomly within a sphere of radius r^^ax — around the centre of V. So 
the probability p(c) for C to be a distance c from the centre of V is 

Mc) = 7^ (14) 

V max) 

Since the normal of this plane has random direction, r can be smaller or equal to c, 
with probability 

p(r;c)f«2^ (15) 

Integrating over all c gives the probability for the plane to be a distance r from the 
centre of V as 

dcp{r; c) p{c) = — — -^{r'^ax - (16) 

Vmax) 

The ensemble average thus becomes an integral over r and s 
MyM = (exp*^0 

r'^ max f ^ X 

= dr p{r) J^^ds -exp{tyi{r)) (17) 
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These integrals can be done, and then eq.(12) can be used to obtain the mgf for X 
which includes the effect of all strings, with the number of strings affecting V at time 
ti being given by 

= KCaJV(^r')' (18) 

There is one problem with the above. The formulas for yi{r) in equations (13) 
are only true if the projection of V onto the plane swept out by the string along its 
normal falls completely into that plane, and is not (partly) outside of it. But the 
latter can happen for large c for some orientations of the plane, since one side of the 
plane, in the direction of the string's motion, has a length 1^ — Hf^ or smaller, so that 
the distance of C to the edge of the plane can be smaller than H~^/2. For li = Hf^ 
one can show that less than half of the strings miss V and give no perturbations to 
it, so that we can estimate this effect by replacing u by i^eff — O.bu . If the strings 
are moving slowly, so that is even smaller, and not at about the speed of light, our 
model is not really applicable because the formula for the imparted initial velocity 
perturbations would change. 

Actually a string can affect V ii c < Ri + H^^ . But for c larger than H^^ we 
encounter the problem mentioned in the previous paragraph, so that we overestimate 
the perturbations by using r^^^x — + ^i- Therefore we calculate the cumulants 
for both r^^j. = Hf^ and r^^^ = + and take their average, and the model 
is more accurate for smaller R. But since at scales below about 5h~^ Mpc nonlinear 
effects become important, and we are only considering linear perturbations, we must 
also keep above that scale. 

4 Probability Distribution and Comparison with 
Observations 

First we want to look at the shape of the probability distribution for Vx- The nongaus- 
sianness is largest on the smallest scales since smoothing makes things more gaussian, 

10 



and larger regions are affected by more strings. Therefore we are going to compare 
p{vx) from strings with the results from IRAS at the smallest scale oi R = 6h~^ Mpc 
considered in Kofman et al. (1994). Also, the derivation of Mx{t) is valid for scales 
of i? < leq, where l^q — 13/i~^ Mpc is the comoving size of the Hubble radius at the 
time of equal matter and radiation. The values Q, — l,h — l/2 and z^q — 2.3 • lO^Q/i^ 
are used. 

Using a symbolic manipulation program (O'Dell 1991), the cumulants are obtained 
from Mx{t) according to cq.(5), giving the following values for the \j (eq.(6)) needed 
in the expansion of the probability distribution (eq.(8)) for R — 6h~^ Mpc : 



These values, as well as the standard deviations quoted below, are the averages of 
the two cases r^^^ = and r^^^ = + i?j . For two values of u, p{vx/cr) is 
plotted in Figure 1, up to the term involving Hs{S) in the expansion (cq.(8)). For 
the higher value of u — 10 strings per Hubble volume, the distribution is practically 
indistinguishable from a gaussian one. For u — 1 there is a slight deviation. We want 
to see if this deviation is significant by performing a -test with the data given in 
Kofman et al. (1994), which has been grouped into bins of size Vx/cr = 0.25. The 
data points fall practically on a gaussian curve, so instead of taking the exact values 
from the data we calculate the absolute frequencies rrij in the j bins expected if the 
underlying distribution were gaussian. The IRAS 1.9Jy survey maps out a sphere of 
radius SOh'^Mpc , so that there are (80/6)^ independent smoothing regions of radius 
6/i^^Mpc. Let rij be the corresponding absolute frequencies expected from the stringy 
distribution. Using 20 inner bins, we find = 4.63 for = 1, where 
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A4 
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0.34/1/ 
0.20/1/^ 



(19) 
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(20) 
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, which is much smaller than the 95% confidence upper hmit of 38.58 for 19 degrees 
of freedom, so that the data is in agreement with the probability distribution from 
strings. If we replace v by z/cfr = v j2 to take into account that the side of the plane 
swept out by the string in the direction of its motion is only half the diameter of the 
Hubble volume, then = 22.5 for = 1. 

Next we want to compare the magnitudes of velocities. The standard deviation 
of the single velocity components is calculated to be 

a = IMv^l'^u for R = 6/i-^Mpc (21) 

(7 = 0.99i/^/^{l for i?=12/i-iMpc (22) 

where u = QAzH^u , and u is defined in eq.(l). We compare these values with results 
from Peacock and Dodds (1994), who used the power spectra of various observations 
to calculate the 3 dimensional rms velocities Vrms of regions of radius R. For a 
gaussian random variable with three independent gaussian variables as components, 
the standard deviation of a single component is given by a = Wrms/v^- Using this 
relation, the values given in Peacock and Dodds (1994) are 

a - (381 ± 156) km/s for R = 6/i"^Mpc (23) 

a = (337 ± 138) km/s for R = 12/i-^Mpc (24) 

where we have taken the fractional error of 1/ y/6 quoted for the actual measurement 
of Vrms at a scale of 5h~^ Mpc. Comparison of a in the string model with the values 
obtained from observations at these two scales, gives as an average value for af 

^ = 3.2 ±0.9 (25) 

where a is the combination of parameters 

" = ^10^^ (26) 
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a can be constrained from the rms value of the temperature fluctuations in the cosmic 
microwave background measured by COBE to be (Perivolaropoulos 1993a) 



q; = 1.0±0.2 



(27) 



Using this value of a , we find 



/ = 3.2 ±1.1 



(28) 



This indicates that there must be some small scale structure on the strings (see 
eq.(l)) in order to obtain the right magnitude of peculiar velocities and consistency 
with CMB observations. This is also in agreement with recent simulations (Bennett 
& Bouchet 1988, Albrecht & Turok 1989, Allen & Shellard 1990), which show the 
presence of small-scale structure on cosmic strings. Our analysis has been done in the 
string wake model, whereas strings with lots of small scale structure accrete matter 
rather in the form of filaments. Therefore the precise value of / is not to be taken 
too seriously. 

Numerical simulations of peculiar velocities from long strings without small scale 
structure have been performed in a similar framework (Hara, Mahonen & Miyoshi 
1993), where u' strings are assumed to move across the horizon at every e-fold expan- 
sion of the universe. The authors found that 



yields good agreement with observations. The number of strings at every two-fold 
expansion used in our analysis is related to i/' by = i/' In 2 . Therefore af — 3.3±0.8 
from these simulations, which agrees quite well with our result of af — 3.2 ± 0.9 . 

5 Discussion 

We have shown that the probability distribution of a single component of peculiar 
velocities in the cosmic string wake model is very close to a normal distribution 



Gil 'JsVs 

10-6 c 



(4±1) 



(29) 
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on scales of several Mpc, as suggested by a general argument for seed models 
(Scherrer 1992), and in agreement with observations. A comparison of the measured 
magnitude of peculiar velocities with that expected from strings, using the normal- 
ization of string parameters from the COBE quadrupole, suggested that strings have 
some small-scale structure. Nongaussian features are more apparent in the velocity 
differences than in the velocities themselves (Catelan & Scherrer 1994), and it would 
be interesting to calculate the probability ditribution of these velocity differences 
within the cosmic string model of structure formation. 
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Figure Captions 



Figure 1: Stringy probability distribution p{S — v^/cr) of a single velocity component 
smoothed over regions of radius R — 6h~^ Mpc for ior u — 1 (solid line) and u — 10 
strings per Hubble volume (dotted line), compared with a normal distribution (dashed 
line) . 
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solid line: nu=1 , dotted line: nu=10, dashed line: Gaussian 
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